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We report the results of a search for a heavy neutral gauge boson Z' decaying into the dielectron 
final state using data corresponding to an integrated luminosity of 5.4 fb _1 collected by the DO 
experiment at the Fermilab Tevatron Collider. No significant excess above the standard model 
prediction is observed in the dielectron invariant-mass spectrum. We set 95% C.L. upper limits on 
a (pp — > Z') x BR(Z' -> ee) depending on the dielectron invariant mass. These cross section limits 
are used to determine lower mass limits for Z 1 bosons in a variety of models. For the sequential 
standard model Z 1 boson a lower mass limit of 1023 GeV is obtained. 

PACS numbers: 13.85 Rm, 14.70 Pw 



The gauge group structure of the standard model 
(SM), SU(3) C ® SU(2) L (g> U(1) Y , can be extended with 
an additional U(l) group, which may arise in models de- 
rived from grand unified theories (GUT) that are based 
on groups with rank larger than four [lj. Additional 
U(l) groups can also arise from higher dimensional con- 
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structions like string compactifications. In many mod- 
els of GUT symmetry breaking, U(l) groups survive at 
relatively low energies, leading to corresponding neutral 
gauge bosons, commonly referred to as Z' bosons [2|. 
Such Z' bosons typically couple to SM fermions via the 
electroweak interaction, and can be observed at hadron 
colliders as narrow resonances through the process qq — > 
Z' — > e + e~ . There is no simple general parametrization 
that can be applied to all the Z' models. Nevertheless, 
the models can be distinguished according to the strength 
of the gauge coupling, gz>, for the additional t/(l) group. 
The models with coupling of electroweak strength are 
called canonical. The sequential standard model (SSM) 
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TABLE I: A selection of commonly used E$ models 

Model sin 8 a 

Z,i, 

z ! 

Z s q 
Z N 



Z' boson is a canonical example, where the SSM Z' bo- 
son (Z' SSM ) is defined to have the same couplings as the 
SM Z boson. The SSM Z' boson is often used as bench- 
mark [2, 3]. An additional example of a canonical model 
can be derived from the superstring inspired Eq mod- 
els Q. The decomposition of Eq can give rise to two 
additional U(l) factors through E 6 -> 50(10) x U(l)^ 
and 50(10) — > SU(5) x U(l) x . These groups are associ- 
ated with the gauge fields Z 1 ^ and Z' x that can mix and, 
at the TeV scale, can give rise to additional Z 1 bosons 
through the linear combination 



Z'{6) = Z'sin9 + ZLcos9, 



(1) 



where < 8 < n is a mixing angle |5|. The most com- 
monly referenced Z' boson models arising from Eq are 
summarized in Table U @ . 

An example of a non-canonical model is the U(l)x 
Stueckelberg extension of the standard model (StSM) 
that gives rise to a very narrow Z' boson 0, The 
Stueckelberg mechanism allows for the possibility of an 
Abelian gauge boson to gain mass without the require- 
ment of a Higgs mechanism. The new parameters that 
are introduced in this model are the StSM mass mixing 
parameter, e, and the Z' boson mass, Mz'- In the limit 
e — > 0, the Stueckelberg sector decouples from the SM 

In this Letter, we report on a search for a Z' boson de- 
caying into an electron pair with the DO detector at the 
Fermilab Tevatron Collider, where protons and antipro- 
tons collide at ^/s = 1.96 TeV. A Z' boson would appear 
as a narrow resonance in the ee invariant mass spectrum, 
with a natural width smaller than the resolution of the DO 
electromagnetic calorimeter. A previous Tevatron search 
by the CDF collaboration [l(|, corresponding to 2.5 fb -1 
of integrated luminosity, sets a lower mass limit on SSM 
Z 1 bosons of 963 GeV and reports a discrepancy over the 
expected SM background at M ee ~ 240 GeV equivalent 
to 2.5 standard deviations. The CDF collaboration has 
also performed a search in the Z' — > /j,fi channel llj . 
corresponding to 2.3 fb^ 1 of integrated luminosity, with 
95% C.L. upper limits on a (pp Z') x BR(Z' /.i^i) 
ranging from ^50 fb to ^3.2 fb for Mz 1 between 175 GeV 
and 1100 GeV. 

The DO detector [13] is composed of a central tracking 
system surrounded by a 2 T superconducting solenoidal 
magnet and a central preshower detector (CPS), a 



calorimeter, and a muon spectrometer. The central 
tracking system includes a silicon microstrip tracker 
(SMT) and a scintillating fiber tracker (CFT) that are 
designed to provide coverage for particles in the pseudo- 
rapidity range \r)\ < 3, where r\ = —In [tan (0/2)], and 9 
is the polar angle with respect to the proton beam direc- 
tion. The azimuthal angle is denoted by 4>. The CPS is 
located between the solenoid and the inner layer of the 
central calorimeter and is formed of approximately one 
radiation length of lead absorber followed by three layers 
of scintillating strips. The calorimeter consists of a cen- 
tral section (CC) covering \r/\ < 1.1 and two end calorime- 
ters (EC) that extend the EM coverage to rj « 4.1, with 
all three sections housed in separate cryostats [l3| . Each 
section consists of an inner electromagnetic (EM) sec- 
tion, and an outer hadronic. The EM calorimeter is seg- 
mented into four longitudinal layers (EMz, i — 1,...,4) 
with transverse segmentation of Arj x Atf> = 0.1 x 0.1, 
except for the finely segmented third layer where it is 
0.05 x 0.05. The muon system, covering |?7| < 2, is lo- 
cated beyond the calorimeter and is composed of a layer 
of tracking detectors and scintillation trigger counters in 
front of 1.8 T iron toroidal magnets, and followed by 
two similar layers after the toroids. The luminosity is 
measured using plastic scintillator arrays in front of the 
end calorimeters. The data acquisition system includes 
a three-level trigger, designed to accommodate the high 
instantaneous luminosity. The data sample was collected 
between July 2002 and June 2009 using triggers requir- 
ing at least two clusters of energy deposits in the EM 
calorimeter and corresponds to an integrated luminosity 
of 5.4 ± 0.3 ftr 
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The event selection requires two isolated electron can- 
didates in the central section of the calorimeter. An 
electron candidate is characterized by an EM clus- 
ter with transverse momentum px > 25 GeV and 
1 77 1 < 1.1, reconstructed in a cone of radius TZ = 

(Ar)f + (A(j)f = 0.4. At least 97% of the EM clus- 
ter energy must be deposited in the EM section of 
the calorimeter and its energy must be isolated in the 
calorimeter, [E iot (0.4) - E BU (0.2)] /E EM (0.2) < 0.07, 
where E tot (TZ) and Eem (TZ) are the total energy and the 
energy in the EM section, respectively, within a cone of 
radius TZ around the electron direction. In addition, the 
EM cluster is required to be consistent with an electron 
shower shape, using a \ 2 test and a neural network dis- 
criminant [15| . The EM cluster is required to be spatially 
matched to either a reconstructed track or a pattern of 
hits in the SMT and CFT consistent with the passage 
of an electron. The scalar sum of the pt of all tracks 
originating from the pp interaction vertex (PV) in an an- 
nulus of 0.05 < TZ < 0.4 around the cluster is required 
to be less than 2.5 GeV. Events are only considered if 
the PV lies within 60 cm of the geometrical center of 
the detector in the coordinate along the beam axis to 
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be fully within the SMT acceptance. The two electron 
candidates are not required to have opposite charges to 
avoid losses due to charge misidentification. The data 
sample consists of 185,264 events that satisfy these se- 
lection criteria in the dielectron invariant mass control 
region 60 < M ee < 150 GeV and 1332 events in the 
search region M ee > 150 GeV. 

Signal and SM background events are generated us- 
ing pythia with the CTEQ6L1 [T3| parametriza- 
tion of the parton distribution functions (PDFs), and 
processed through the DO detector simulation based on 
GEANT3 [HI adding zero bias events, and the same recon- 
struction software as the data. Signal templates based on 
the SSM Z' boson have been generated up to masses of 
1100 GeV. The width of the resonance scales with the Z' 
boson mass, according to Tz> = Tz x Mz'/Mz, where 
Mz and Tz are the mass and width of the Z boson. For 
Mz> > 2mt the decay channel to top quarks opens up, 
thus increasing the width of the resonance. The signal se- 
lection efficiency increases from ~22% to ~44% for Mz< 
between 175 and 1100 GeV independent of the type of 
Z' boson discussed in this Letter. 

The dominant irreducible background is due to the 
Drell-Yan (DY) process. A mass-dependent k-factor [l9[ 
has been applied to the pythia dielectron invariant mass 
spectrum to account for next-to-next-to-leading order 
(NNLO) contributions. The main instrumental back- 
ground originates from the misidentification of one or 
two jets as electrons. The shape of the invariant mass 
spectrum for this background is obtained from data by 
selecting events where the EM clusters fail the x 2 test- 
Other SM backgrounds include Z/7* -> tt, W+j, WW, 
ZZ, WZ, W+ jets, tt, and 77 production. The contri- 
bution of these background processes is small (~0.6%) 
and is estimated using PYTHIA corrected for higher order 
contributions ^d^. 

The normalization of the various background con- 
tributions is determined by fitting the invariant mass 
spectrum of the data to a superposition of the back- 
grounds in a control region around the Z boson mass 
(60 < Mee < 150 GeV), where the existence of Z' bosons 
has been excluded by previous searches 0. The to- 
tal number of background events in that region is fixed 
to the number of events that have been observed in the 
data. The relative contribution from the DY process and 
instrumental background is a free parameter, while the 
contribution from the other SM processes is normalized 
to their theoretical cross sections. The uncertainty of the 
background normalization is estimated by varying both 
the criteria to select the instrumental background sample 
and the fitting range, and is 2%. 

Having normalized the various background contribu- 
tions to data in the control region, the background shapes 
are used to extrapolate to higher invariant masses. The 
measured ee invariant mass spectrum, superimposed on 
the expected backgrounds for the full mass range studied, 
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FIG. 1: Distribution of M ee for data, along with the total 
expected background for the full invariant mass range studied. 
Insert focuses on the area of the M ee spectrum from 160 GeV 
to 300 GeV, where the majority of observed data in the signal 
region lie. 

is shown in Fig. [TJ The data and expected background 
are generally in good agreement for the full invariant 
mass range studied, with a y 2 over degrees of freedom 
equal to 118.5/113. 

In the absence of a heavy resonance signal, the ee in- 
variant mass distribution is used to calculate an upper 
limit on the production cross section of Z' bosons multi- 
plied by the branching ratio into the ee final state, using 
a Poisson log-likelihood ratio (LLR) test statistics [24j . 
The expected limits are calculated using the median of 
the LLR distribution for a background-only hypothesis. 
The observed limit, obtained including all the fluctua- 
tions present in the data, is expected to be contained in 
the ±1 and ±2 standard deviations region with a prob- 
ability of 68% and 95%, respectively. An observed limit 
significantly outside the expected range would indicate 
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FIG. 2: The observed and expected 95% C.L. upper limits on 
a (pp — > Z') x BR (Z' — > ee) as a function of M z > , compared 
to the theoretical predictions of the cross section for the SSM 
Z' boson and the Z' bosons arising from the E§ model. The 
median expected limits are shown together with the ±1 and 
±2 standard deviation bands. 



either a poor modeling of the background or that the 
data is inconsistent with the background-only hypothe- 
sis. 

The following systematic uncertainties on the expected 
background and the signal have been considered for the 
limit calculation. The uncertainties affecting the ex- 
pected background include the electron identification ef- 
ficiency (3.0% per electron), the mass dependence of the 
DY associated NNLO k-factor (5.0%), and the back- 
ground normalization (2.0%). Uncertainties that affect 
the signal include the integrated luminosity (6.1%), the 
PDFs for signal acceptance (0.4% - 7.6%), the electron 
identification efficiency (3.0% per electron), the EM clus- 
ter energy resolution (6.0%), and the trigger efficiency 
(0.1%). For the EM energy resolution and the back- 
ground normalization, both the effects on the normal- 
ization and on the shape of the invariant mass distri- 
bution have been considered in extracting limits. For 
the remaining systematic sources only the changes to the 
overall background normalization or signal detection effi- 
ciency have been considered. The systematic uncertain- 
ties are included via convolution of the Poisson probabil- 
ity distributions for signal and background with Gaus- 
sian distributions corresponding to the different sources 
of systematic uncertainties taking into account all rele- 
vant correlations between systematics' sources. 

The observed upper limits on the production cross sec- 
tion multiplied by the branching ratio into an ee pair for 
the process pp —> Z' — >• ee are given in Table [TTI as a func- 
tion of the mass hypothesis, together with the median ex- 
pected limits calculated under the assumption that the 



FIG. 3: Excluded region in the (M z ,, g z >) plane at 95% C.L. 
for the Z' x model. The expected limit is superimposed. 



observed dielectron invariant mass spectrum arises only 
from the backgrounds considered in the analysis. Fig- 
ure [2] shows these limits together with the ±1 and ±2 
standard deviation bands on the expected limit, and the 
cross section predictions for SSM and Eg Z 1 bosons 0] 
where a constant k-factor of 1.3 25] has been applied to 
the PYTHIA cross section. Since this analysis searches for 
a resonance instead of an enhancement in the total cross 
section, signal cross section predictions are calculated by 
integrating over the region [Mz> — 10 x Tz>, 00], where 
Tz> is the width of the SSM Z 1 boson, thus excluding 
Z' boson events which do not contribute to the resonant 
region. For Mz' < 500 GeV the difference between the 
cross section in the region denned above and the total 
cross section is less than 5%, while for a Mz' = 1 TeV 
SSM Z' boson it is ^40%. The mass limits on the specific 
models of Z' bosons considered are given in Table Mil 

These limits can be translated into upper limits on the 
U(l)z> gauge coupling, g z > @, as a function of Mz 1 - Fig- 
ure [3] illustrates the observed upper limits on gz 1 / gz' [26| 
for the Z' x model. 



Cross sections are calculated as a function of Z' bo- 
son mass to interpret the observed upper limits on 
a (pp -> Z') x BR (Z' -> ee) as mass limits for a StSM Z 1 
boson. Figure 0] shows the observed and expected limits 
and the cross section predictions for the StSM Z' boson 
for several e values from 0.02 to 0.06 @. The mass limits 
are summarized in Table [TlTl 

In summary, we have searched for a heavy narrow res- 
onance in the ee invariant mass spectra, using 5.4 fb _1 
of integrated luminosity collected with the DO detector 
at the Fermilab Tevatron Collider. The observed spec- 
trum agrees with the total background expected from SM 
processes and instrumental backgrounds. No evidence 
for physics beyond the SM is observed. For a Z' boson 
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TABLE II: Expected and observed 95% C.L. upper limits on 
the production cross section multiplied by the branching ratio, 
a (pp -> Z') x BR (Z' -> ee). 

(fb) 



TABLE III: Expected and observed lower mass limits for var- 
ious Z 1 bosons. 



Z' Boson Mass 


rr ( tin — ^ 
u \PP ' 


Z') x RR (Z' — > ee) 


(GeV) 


Expected 


Observed 


175 


49 


22 


900 


o u 




99^ 

Lid 


9Q 


1 6 


250 


23 


15 


97^ 


9(1 


1 


300 


1 6 


24 






90 






7 


375 


10 


6.9 


400 


o . o 


7 9 


450 


fi 8 


8.2 


500 


5.5 


5.4 


550 


4.4 


6.2 


600 


3 7 




DOU 


1 


•-t 


700 


2.7 


3.2 


750 


2.4 


3.2 


800 


2.2 


2.6 


850 


2.2 


2.3 


900 


2.0 


2.1 


950 


1.9 


2.0 


1000 


1.9 


2.0 


1050 


1.9 


1.9 


1100 


1.9 


1.9 



D0, 5.4 fb 1 




q<w ri J ~ Observed 

9b/oU -l Expected 

Expected ±1a 
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— SSM Z' 

— Z stSM e = 0.02 

Z StSM 

^ StSM e 
— Z ,StSM e 
StSM * 



0.03 
0.04 
0.05 



200 300 400 500 600 700 800 900 10001100 

M z , (GeV) 

FIG. 4: The observed and expected 95% C.L. upper limits on 
a (pp — > Z') x BR (Z' — > ee) as a function of Mz> , compared 
to the theoretical predictions for the Z' boson cross sections 
in the SSM and in the StSM extension for values of e ranging 
from 0.02 to 0.06. The median expected limits are shown 
together with the ±1 and ±2 standard deviation bands. 



Model 


Lower Mass Limit (GeV) 




Expected 


Observed 


~7 

^SSM 


1024 


1023 


Z n 


927 


923 


Z \ 


910 


903 


Z j, 




on i 


Z^ 


879 


874 


Z'sq 


829 


822 


Z'x 


795 


772 


ZstSM(e = 0.06) 


471 


443 


^stSM(e = 0.05) 


414 


417 


Z'stsmit = 0.04) 


340 


289 


Z' B tSu(. e = 0.03) 


227 


264 


Z' StSM (e = 0.02) 




180 



with SM couplings and with intrinsic width significantly 
smaller than the detector resolution, we set 95% C.L. up- 
per limits on a (pp — > Z') x BR (Z' —> ee) between 22 fb 
and 1.9 fb for M z , between 175 GeV and 1100 GeV. 
These represent the most stringent constraints to date, 
and translate into a lower limit on the mass of the SSM 
Z' boson of 1023 GeV. 
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